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NOTATION 
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The  concentration  of  air  in  the  water  at  an  arbitrary  po.n- 
and  time 

The  concentration  of  air  in  the  water  at  infinity;  also  the 
concentration  of  air  in  water  when  the  water  is  in  equilibrium 
with  air  at  pressure  P0  and  temperature  T0 

The  specific  heat  for  air  under  constant  pressure 

The  specific  heat  for  air  under  constant  volume 

The  diffusion  constant  for  air  in  water 

f = ^«-v^-1^coeV^-l) 

G l) 

/ (V)  =*  f e~v^‘  ^ ar~*sm  V(ar— l)d* 

0 J i 

The  Boltzmann  constant 

Tt.e  kinetic  energy  of  the  motion  of  tr.e  water 
The  average  molecular  mass  for  air 
Number  of  molecules  in  tr.e  bubble 
The  initial  number  of  molecules  in  ti.c  oubble 
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Tr.e  rate  at  wr.icn  molecules  leave  tr.e  bubole  and  enter  the 
water 

Tr.e  rate  at  wn’.ch  molecules  stream  away  from  water  sur- 
face anc  enter  tree  buhcle 

The  molecular  density  of  the  air 

Tr.e  molecular  density  for  air  at  pressure  P0  and  temperature 

r. 

The  pressure  in  the  bubble 

The  undisturbed  hydrostatic  pressure;  the  average  pressure  at 
infinity 

The  sound  pressure 

The  pressure  in  the  water  at  infinity 

The  relative  amplitude  of  tr.e  sound  pressure;  p sin  ut 

The  instantaneous  radius  of  tne  bubble 

The  time  rate  of  change  of  R 

The  equilibrium  radius  of  the  bubble 

The  resonant  radius  at  angular  frequency  to 

The  distance  from  the  center  of  the  bubble  to  a point  In  the 
water 

The  temperature  of  the  air  in  bubble,  Kelvin 
The  constant  temperature  of  the  water,  Kelvin 
Time 

The  radial  velocity  of  the  fluid  at  the  point  r 


The  rate  at  which  work  is  dene  upon  the  water 


I 


- (m/ZkTV* 

The  value  of  0 at  T - Tc 

* The  ratio  of  specific  heats  • approximately  4/3  for  air 

e, 

at  standard  conditions 
Damping  constant 
A constant  of  integration 


A = (-A-V  \ & 
\2irpJ  Ca.  D 


The  density  of  the  water 

The  density  of  air  at  pressure  /*  and  temperature  T% 

The  period  of  the  sound  beam 
Frequency 

Natural  frequency 
Angular  frequency 

The  angular  frequency  corresponding  to  natural  frequency 


l 


THE  DEAERATION  OF  WATER  BY  A FOUND  BEAM* 


vz  rode,  X.F. 


ABSTRACT 


A calculation  is  made  to  determine  what  part  might  be  played  by  the  diffusion  of 
air  into  existing  bubbles  as  a consequence  of  the  mechanical  motions  Induced  fcy  a sound 
beam.  The  case  of  an  isolated  bubble  iu  a weak  sound  beam  of  wav<;  length  considerably 
greater  than  the  bubble  radius  is  considered  Th“  effects  of  surface  tension,  viscosity  and 
energy  dissipation  are  neglected.  The  diffusion  problem  Is  treated  by  means  rf  a per- 
turbation technique.  The  growth  of  the  bubble  is  found  to  be  of  second  order  in  the  sound 
pressure, l.e.  proportional  to  the  sound  intensity.  Numerical  results  presented  show  that 
the  effect  is  sufficient  to  account  for  significant  bubble  growth  especially  in  the  case  of 
very  small- bubbles. 


INTRODUCTION 

The  deaeration  of  water  by  an  ultrasonic  beam  Is  a phenomenon  which 
has  become  Increasingly  familiar  to  experimenters  In  ultrasonics.  However, 
there  Is  as  yet  no  satisfactory  account  of  the  mechanism  responsible  for  this 
effect.**  The  calculation  that  is  made  here  lends  support  to  the  view  that 
the  deaeration  results  from  the  growth  of  small  air  bubbles  which  Is  caused 
by  tne  mechanical  motion,  that  Is  forced  oy  the  souna  beam.  The  motion  of  a 
bubble  responding  to  a sound  beam  Is  unsymmetricai  In  suen  a mannei  that  trie 
average  raclus  becomes  slightly  greater  tnan  the  equilibrium  radius.  More- 
over, because  tr:e  diffusion  of  a'r  ' nto  Me  bubble  varies  with  the  area  of 
the  boundary,  the  periods  when  the  bubble  Is  expandea  have  greater  effect 
upon  the  ciffusion  of  air  Into  the  bubble  -han  the  periods  when  the  buoole 
is  contracted,  hence  * t is  to  be  expected  that  the  mechanical  motions  of  an 


•A  dissertation  suhtitter1  to  the  f: tcr-lty  ■ : thx  Gscduatc  School  of  Arts  and  Sciences  of  the  Catholic 
University  of  America  iu  parti*)  fulf : lit cut  of  the  requirements  for  the  degress  of  Doctor  of  Philosophy. 

••After  most  of  the  work  prevented  here  read  seen  accomplished  it  was  discovered  that  a similar  cal- 
culation had  b-en  made  by  F.C.  Plaice.  Ur.  cie  fefs-er.ee  2 on  page  20  ol  this  report.  A comparison 
of  the  two  calculations  )s  presetted  lr.  Appendix  I of  rr.es  paper. 
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air  bubble  will  result  in  the  growth  of  the  bubble  Tr is  mechanism  for  the 
growth  of  an  air  bubble  in  a sound  field  is  very  closely  related  to  the  mecn- 
i sm  for  growtr.  of  an  air  bubble  in  a cavitating  i-.ow.  j.ic  xatter  case  has 
been  treated  in  an  earlier  paper.1  The  main  difference  between  t.'.ese  two 
phenomena  are  the  much  smaller  scale  and  greater  symmetry  oi  tr.e  motion  of  t:ie 
sound  field.  On  the  other  hand  a bubble  can  readily  be  exposed  to  a sound 
f'eld  for  a relatively  long  time  whereas  tr e time  spent  in  a cavitating 
flow  Is  usually  brief.  The  calculation  that  is  presented  here  follows  along 
the  lines  of  the  earlier  calculation.  However,  because  of  the  small  scaie  of 
the  motion  it  is  now  possible  to  treat  the  diffusion  of  air  threngr.  the  water 
in  more  detail . 

A single  small  bubble  is  considered  to  be  immersed  in  an  infinite 
body  of  water  through  which  a sound  beam  is  passing.  The  rad* us  of  the  bubble 
is  assumed  to  be  so  small  in  corapar'son  to  the  wave  length  of  the  sound  beam 
that  the  beam  can  be  represented  by  a periodic  variation  ir.  the  r.ydrostatic 
pressure  at  infinity.  Consequently  spherical  symmetry  is  preserved  'n  the 
motion.  Only  the  case  of  a weak  sound  beam  is  considered  so  that  the  varia- 
tion of  the  pressure  at  infinity  may  be  considered  small  compared  to  the  aver- 
age hydrostatic  pressure.  Surface  tension  is  neglected  and  it  is  assumed  that 
all  conditions  are  such  that  in  tne  absence  of  the  sound  oeam  equ* 1 !br* urn 
would  obtain.  The  conduction  of  neat  is  also  neglected  so  tnat  tne  tempera- 
ture in  the  water  remains  constant  and  uniform  ar.a  the  compressions  ana  ex- 
pansions of  the  air  in  the  bubble  are  adiabatic.  In  addition  the  preser.ee  of 
water  vapor  in  the  bubble  is  ignored. 

The  ooject  of  this  calculation  is  to  determine  the  average  rate  of 
growth  of  the  bubble  wnen  conditions  have  become  steady.  It  will  be  assumed 
-hat  the  change  of  the  amount  of  air  in  the  bubble  during  t •■.mu  required 
to  attain  steady  oundit'ons  ana  during  a period  of  the  sound  fie. a are  botr. 
insignificant  in  comparison  to  the  amount  of  a*r  initially  the  bubble. 
Therefore,  the  effect  of  the  change  of  air  content  during  one  period  on  the 
rate  of  growth  will  be  ignored.  Thus  the  pressure  sf  the  -'i-  bubble 

will  be  assumea  to  be  directly  given  by  tne  adiabatic  relations.-. ‘p  w i : nout 
correst'on  for  the  change  in  the  air  content  in  the  buboie.  J'miiar.y  t;.e 
density  of  air  In  the  bubble  will  oe  taken  as  inversely  presort* cnal  tc  the 
volume. 

The  calculation  is  d'videb  into  three  parts  covering  (1)  s,.e  equation 
of  motion,  of  the  boundary  of  the  bubble.  [2)  tne  boundary  conditions  for  the 
diffusion  equation,  C'.j  tne  solution  of  the  d f'u-*  • 

^fl«for»nco»  &r«  on  pa^e  20. 
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THE  EQUATION  CF  MOTION  OF  THE  BOUNDARY  OF  THE  BUBBLE 

The  equation  of  motion  of  an  air  bubble  in  water  under  the  conditons 
spec*fied  above  has  been  treated  by  many  writers.  The  derivation  is  repeated 
nert  in  order  to  obtain  the  equation  in  a form  suitable  for  present  purposes. 

Let  the  origin  of  coordinates  be  the  center  of  the  bubble.  Let-  r 
be  the  distance  from  the  origin  to  a point  in' the  water  and  let  t>  be  the  radi- 
al velocity  of  the  fiuid  at  this  point.  Let  R be  the  radius  cf  the  bubble  and 
ft  be  the  time  rate  of  change  of  R.  Then  the  incompressibility  of  the  water 
requires  r*  . 

v=7*R  [l] 

and  the  rate  of  change  of  kinetic  energy  is 

y4*r*ir- « />.[*«*■ +!*’*’]  m 

where  Pw  is  the  density  of  the  water.  If  Pt  is  the  pressure  in  the  bubble  and 
JJ,  is  the  pressure  in  the  water  at  infinity,  the  rate  at  which  work  is  done 
upon  the  water  is  given  by 

(r.-vV-ir  i?> 

Writing  Pv  = PQ  4-  Pt  where  Pt  is  the  undisturbed  hydrostatic  pressure  and  P, 

1 s the  sound  pressure,  and,  since  the  motion  is  assumed  to  be  adiabatic, 
p.”  M . y being  the  ratio  of  specific  heats,  Cp/e,  ; for  the  air  in 

the  bubble,  this  equation  becomes: 


Since  no  damping  effect  is  included  dKE/it  - dW/dt  and  the  general  equation  of 
motion  may  be  written 


**  +4 ». «’ + p,  fi  - ( f*f]  + p.-» 


For  sonic  motions  we  may  write  R = /?0+  A R , where  is  the  equi- 
librium radius  of  the  bubble  and  AR  <<  ftt  whence  to  second  order  in  AR  the 
equation  of  motion  [5]  becomes: 

A./?DAft  + 8rp0-^+ftp4/?i/?+|/uU«),--!y(3y  + i)p*(^),+  p*“0  M 

To  first  order  for  free  sonic  oscillation, i.e.  P,  * 0 the  familiar  result  is 
obtained 

2 , 2 2 3yp0 

w02=  v0  - — pr  r 71 


[73 


where  v0-^is  the  natural  frequency  and  u)0  the  corresponding  angular  frequency. 

Writing  P,  = p P0 sin u t f where  p<<l  the  following  solution,  to  the 
second  order  in  p,is  obtained  by  the  standard  perturbation  technique.  This 
solution  can  be  directly  verifiea  by  substitution  in  Equation  (c) 


where 


AR  — R — Rq  = rq  a0(sin  ut)p+(a,  + a2eos  2 uit)  p2 


/ cj2  \ 


3>’  + 1 


37  + 1 _ " 

a2(  ) * A 

0 ' 4 I 64 


3y  + l + 5-^ 



3Y  + 1 + 5 


4(3+(l-!U)  (1-4^)  4(l-4^) 


The  numerical  approximations  are  obtained  by  using  y - 4/3  for  air  and  ignor- 
ing (j/u0  when  compared  to  unity  in  accordance  with  the  hypothesis  that  we  are 
dealing  with  a bubble  whose  radius  is  small  compared  to  the  radius  for  reso- 


Differentiation  of  Equation  [8j  gives 

R — R0  ( a0  up  cos  u)  t - 2a2a>p2  sin  2 


THE  BOUNDARY 'CONDITION  FOR  THE  DIFFUSION  EQUATION 

The  traffic  of  air  molecules  at  the  boundary  of  the  bubble  is~studl- 
ed  by  considering  the  rate  dNJdt  at  which  molecules  leave  the  bubble  and  enter 
the  water  and  then  the  rate  d Afc/di  at  which  molecules  stream  away  from  the 
water  surface  and  enter  the  bubble  thus  obtaining  the  net  rate  of  transfer  of 
molecules  ~ dNJii  which  is  the  rate  of  growth  of  the  bubble 

Because  of  the  spherical  symmetry  it  is  only  necessary  to  find  j+p  rate  of 
transfer  of  molecules  per  unit  area  of  the  boundary  (rfAfc  f)  + 4***  This  rate 
may  ce  determined  by  considering  just  a small  section  of  the  boundary. 


,v”.<  i-i 


It  is  assumed  that  the  mean  free  path  of  the  molecules  in  the  bubble 
is  sufficiently  small  compared  to  the  radius  of  the  bubble  so  that  the  section 
of  the  spherical  surface  may  be  taken  small  enough  to  be  virtually  a plane. 
From  kinetic  theory  the  rate  of  flow  of  molecules  per  unit  area  through  a 
plane  surface  element  in  the  bubble  is  given  by  n/2K5ryS  wnere  n is  the  mole- 
cular dens' ty  and  ~ Ym/2kT  ,m  being  the  molecular  mass,  T the  temperature 
(Kelvin)  of  the  air  in  the  bubble  and  k the  Boltzmann  constant.3  Since  the 
motion  is  assumed  to  be  adiabatic  /3~ ^(R/r^2^-1)  where  j80  Is  the  value  of 
/?  at  equilibrium  i.e.  A>  = (m/2Jc^)V2=-  i^/ZPo)^2,  and  pa  is  the  density  of  the  air 
equilibrium.  Also  «=  n0(R0/R?  where  Is  the  molecular  density  at  equi- 
librium. Although  some  of  the  molecules  that  strike  the  surface  of  the  water 
rebound  into  the  gas,  almost  all  remain  in  the  water  and  It  is  sufficiently 
accurate  for  our  purpose  to  assume  that  all  the  molecules  that  strike  the 
boundary  remain  in  the  water.*  Thus  the  specific  rate  of  loss  is 


_ „ 1 n V 

°\2np„ I \ R I 


— r.  — (V  + l) 

2 (RqV(  ’ 


The  rate  at  which  molecules  stream  into  the  bubble  from  the  water 
surface  depends  upon  the  concentration  of  the  dissolved  air  and  the  temperature 
of  the  water.  The  density  of  air  in  equilibrium  with  water  when  both  are  at 
the  same  fixed  temperature  Is  proportional  to  the  concentration  of  dissolved 
gas  in  accordance  with  Henry's  Law.  The  effect  of  temperature  is  to  vary  the 
constant  of  proportionality.  Let  c(r,t)  be  the  molecular  concentration  of 
air  In  water.  Then  cx  = c(°°,t)  is  the  molecular  concentration  of  air  In 
water  when  the  water  is  in  equilibrium  with  air  at  pressure  P0  and  temperature 
T0  . Since  the  temperature  of  the  water  is  everywhere  fixed  at  T0  , when  the 
concentration  of  air  in  the  water  at  the  boundary  of  the  bubble  is  cB  the 
rate  at  which  molecules  leave  a unit  area  of  the  water  surface  is  no(Io/2iTft1)1^' 
When  the  concentration  of  air  in  the  water  at  the  boundary  Is  c(R,t ) the 
water  would  be  in  equilibrium  with  air  at  temperature  T0  and  molecular  density 
(*(* .t)/cK)n0  so  that  the  specific  rate  at  which  molecules  stream  away  from  a 


♦This  may  not  be  the  case  when  impurities  have  collected  at  the  boundary  of  the  bubble.  It  has  been 
observed  for  example  that  the  deaeration  of  tap  water  that  has  been  quiescent  for  a long  period  of  time 
will  not  occur  until  the  sonic  pressure  is  increased  beyond  a minimum  threshold  level.  The  suggestion 
has  been  made  by  F.  E.  Fox  and  K.  F.  Herzfeld  that  the  accumulation  of  an  organic  "skin"  at  the  boundary 
of  the  bubble  may  be  responsible  for  this  phenomena.  Such  a skin  would  stifle  the  diffusion  process. 

The  calculation  made  here  presupposes  a fresh  surface  at  bubble  or  water  free  of  impurities  with  the 
f»shold  pressure  is  obtained. 
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unit  area  of  the  water  sui  'ace 


= , ( JLj 

4'lrr  R*  °V2t rpj 


P0  f cm) 
2n  pj 


The  net  rate  ol  transfer  is  thus  given  by 


The  net  rate  of  transfer  ol  mo'.eculis  through  the  boundary  i>  alsv 


V V f y = ■-  4ttR'd(~) 

di  \ d T ' r _ 


wiit.'e  D is  the  diffuse  or  ednstar.t  for  air  in.  water  and  c is  giver; - 'n  ';■■■ 
molecular  concentration'  Esnes  "he  boundary  ;ondiliou  at  r R dt  be 


fl,L  ( is ) A f«  .R»  <>  _ ! 

c o \ l>t  !,■*/"  L too  R> 


whe'-e 


^ ha?  mm 


Under  normal  conditfcas  A in  exceedingly  large  For  standard  cone \t ’on;  > 
temperature  and  pressure  thb  following  values  nuy  be  used; 

•'  P | 5 i ■ ''  : d . : i \ ■'  ;V 

P0  = 1.01  x 106  dynt  cm-2  - D = 1.9  x 10" 6 cm2  f?c_1  ; Cce  = 5.02  * 1017  ?w  ~3 

Wo  =■-  2.50Jrl()19  cm"3  ; pa=  1.20x  J0“3  gm  cm'3. 

Then  for  a bubble  who  S3  equ! librium  radius  Is  lb  3 cm  , A ^ 3 x 1 06 . 

For  all  cases  of  ore, sent  interest  the  value  of  the  left  hand  side  of  hq  .a1  '01 
[ 17 i is  always  small!  relative  to  A.  Hence  little  error  is  made  if  the  equ  IHt 
is  divided  by  A and  set  equt 1 to  zero  so  tna/.;  the  boundary  condition  bee -me  . 

„ fD  >\ ■ . 1 —5 

C l <■«  ( [,  l ' . f If  1 

’•it  r 

This  condition  is  s'-. -d y^  the  statist. it  r'.ia^  ,e  '..icei,  .•A.'i.y.'  jissolvec  air 
in  the  water  at  the  bound..  ~-af_the  bob!  -e.  instar-'^^islv  5s sum*-,  one  value 
required  for  equilibrium  wr:h  ^ ...  . ......  .,/•  t.  ; r, the  boundary 

condition  is  used  in  the  form  of  Equation  h ie j 1;  «E  essential  that  the  rate 
of  growth  of  the  bubble  be  computed  with  the  formula  given  by  Equation  [l6j 
and  not  with  that  of  Equation  Ll5j  since  the  latter  would  Involve  the  product 
of  a very  large  number  and  s very  small  number  whi  :h  has  been  set  equal  to 


7 


zero  in  Equation  [18], 

THE  SOLUTION  OF  THE  DIFFUSION  EQUATION 

The  diffusion  of  the  dissolved  air  through  the  water  is  described 
by  the  partial  differential  equation: 


. „ , •>  . D d / , dc  ^ R~  6 

(DVZ  - v-v)e  — r,gVr  I- 


dc  _ d c 
dr  d t 


. here  and  R are  given  as  functions  of  time  by  Equations  [8]  and  (12]  re- 
-i;  ectiuely . The  boundary  conditions  are  Imposed  by  the  fixed  concentration 
at  infinity,  i.e.  c(°°,l)=ca>  and  either  Equation  [1?]  or  Equation  [ 18  ] . 
Although  there  is  no  inherent  difficulty  in  employing  the  more  precise  Equation 
t 17] , the  labor  involved  in  using  Equation  [l8]  is  considerably  less  and  in 
view  of  the  approximate  nature  of  this  calculation  in  other  respects,  little 
i:>  lost  by  the  use  of  the  Equation  [l8j.  In  either  form  the  boundary  condition 
al  r = ii  is  in  general  awkward  to  handle  because  the  condition  is  imposed  at  a 
moving  boundary. * However,  this  difficulty  can  be  avoided  In  the  present  case 
because  the  motion  of  the  boundary  is  small  so  that  equivalent  conditions  at  a 
flsed  boundary  i.e.  r = can  be  extrapolated  from  the  condition  at  the  moving 
boundary . 

A solution  is  sought  which  is  a perturbation  from  the  undisturbed 
stale.  Thus  c(r,t)  is  expanded  In  powers  of  p : 

c(r,t)=  cx  + cHr,t)p+  c2(r,t)p2  + - - [20] 


and  st r = R expansions  of  the  following  type  are  employed:  (see  Appendix  II). 
c(.im)  = c(R0,t)+(R-R0)~(RQit)  + j(R~R0)2  —2  (R0,t)  + ~ 121] 


/#In  c:)nilt  'tion  with  a similar  heat  diffusion  problem,  M.  S.  Plesset  and  S.  A.  Zvick?,  have 
~--^%c.loyed  £ t schnique  for  treating  the  condition  at  the  moving  boundary  that  could  be  adapted  to 
' v preseat  'roblem.  This  technique  is  based  upon  the  fact  that  the  diffusion  takes  piece  mainly 
^ the  wa  .61  lTanediately  surrounding  the  buoble.  Successive  approximations  can  be  made  each  ef- 
- lively  ei  ending  the  diffusion  to  a larger  volume  of  the  fluid.  However  in  the  present  case 


~.w  . — — — — — ~ ^ rir  wwiaj.u. 

b®  P Jdioated  upon  the  smallness  of  the  motion.  The  method  employed  In  this  paper  U3es  tlie 
that  or  . small  motions  are  to  he  considered  to  avoid  any  approximation  in  regard  to  the  region 


'..;--..^3in  the  1 f fusion  takes  place. 


• irra»-  - r — ^a1*1  «r*wrv<y,'E  um* 


[2>]  in  [18].  Thu: 


e^J)_dc^)  +(R_/?o)p(^(<)+|(R-R0)J^(/?0,0+  ---  [22] 

The  assumption  that  equilibrium  exists  In  the  absence  of  the  sound  beam  re- 
quires that  the  first  terra  in  the  expansion  of  c in  Equation  [ 20 J be  iden- 
tically Cm  as  written  and  it  Is  easily  verified  that  c(r,t)  = c„ 
satlfles  the  differential  equation  and  the  boundary  conditions  with  the  ex- 
clusion of  terms  of  and  beyond  first  order  In  p . Using  Equations  [8]  and  [20] 
in  Equations  [21]  and  [22]  one  obtains  to  second  order  in  P 

c{R,t)=  c„.+  Cl(R,t)p  + [R0o^( sin  <Jt)  («,,()+  <?,(*<,, <)]p*  t23] 


0 C { R § t)  d C i (Rf),t)  , f n /•  j \ ^ ^1  / n i\  . ® ^2  /n  .y1*  l 'll.  ' 

-Q- -g~-^  - P +[R0a0(sinwO^(R0.<)+  ^T^.Ojp*  [2^j 

Substitution  from  Equations  [8]  and  [24]  into  Equation  [l6]  and  the 
taking  of  time  averages  results  in  the  following  expression  for  the  average 
rate  of  growth  of  the  bubble  to  second  order  in  p 


[23] 


and  the  boundary  condition  at  infinity  is  simply  <:,(<»,£)-  0.  Hence  the  de- 
termination of  c,  is  mathematically  identical  to  the  one-dimensicnai  problem 
cf  heat  conduction  in  a semi-infinite  medium  with  prescribed  time-dependence 
of  the  temperature  at  the  boundary.  The  solution  of  this  heat  conduction  pro- 
lem  is  well-known.4  Application  of  this  known  solution  gives 


where 


and 


C1(J  * *L  Oo  [Fan  ut  - G cog  tat] 

F « e~  v (*o  ~ ‘)  cos  v( £ - l) 
r '«o  ' 

G =^e-K«o~1)8inV'(-f-l) 
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This  solution  can  also  be  readily  obtained  by  expansion  of  c,(r,t)  In  a 
Fourier-serles  in  t . 

The  following  evaluations  can  now  be  made: 


jj£p+  2a0(ain  wr)^J  + ^a0(sin  w t)  l |l2j  p 2 


[2}] 


where  all  derivatives  with  respect  to  r are  to  be  evaluated  at  r = R0  . 

It  is  apparent  from  Equation  [2}J  that  although  the  time  dependence  of  Cj  is 
required  for  the  determination  of  the  average  rate  of  growth  to  second  order 
in  p only  the  time  average  of  c,  is  needed. 

The  differential  equation  for  c,  is  found  by  applying  Equations  [8], 
[12]  and  [20]  to  Equation  [19],  expanding  in  powers  of  p and  equating  the 
coefficients  of  the  first  power  of  p . Thus 


d d / 2 ac,\  _ a_c 

r2  §r'  fir'  d t 


dc,\  d c, 


or  the  equivalent 

l 


D§T2 


[26] 


[27] 


The  boundary  condition  at  r - R,  is  likewise  obtained  using  Equations  [8]  and 


fic,  A ,0c,  3V  + 3 (1  + V)  . 9*c. 

i.  ~ rt  . Qin  M t ■ s — n y*  _ __l  • am  ^ — — V — i 


= 0 ; sin  w t 
dr  Or 


4 aoc’ 


; sin 


Or* 


The  differential  equation  for  is  obtained  by  applying  Equations 
[8],  [12]  and  1 20 J to  Equation  [19].  expanding  in  powers  cf  p , equating  the 
coefficients  of  the  second  power  of  p, substituting  the  result  Just  obtained 
for  c,  and  taking  time  averages.  Thus  one  obtains 
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The  boundary  conditions  at  r is  likewise  obtained  using  Equations  [8] 

and  [23]  in  [l8]. 


* According  to  tho  Interpretation  of  diffusion  as  a conaequencs  of  the  randoe  wandering  of  aolecules 
as  presented  by  Klnateln,  ve  nay  write  D *=*V2 r where  **  ie  the  average  square  of  the  displacement 
of  an  air  nolecule  in  vater  after  the  tiae  r which  aay  be  taken  as  the  period  of  the  sound  beea.  Then, 
since  «»  2 iffy  , V m = (2  #o  that  V nay  be  interpreted  a iruitjpie  of  the 

ratio  of  R,  to  the  root -nean- square  displacement  in  a period  of  the  sound  beaa. 


'■•’PW: 
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L - . r — ■ '~*V-  -='  • - " '•'■ - i:'"?" it*** 


5i^)  = .(»t±l)[(£±i-^)=,*  + 0l]  i«i 

Equation  [30]  Is  easily  Integrated  to  obtain 

+ ^ I32J 

where  * Is  a constant  of  integration ; .a  secona  constant  of  integration  vnn1sr.es 
because  of  the  boundary  condition  at  infinity,  c2(<»)  s o • Substitution 
from  Equation  [32]  into  Equation  [31]  permits  the  evaluation  of  *■ 

‘ 2 / 0 L8  8 + 2 Aojjj  V \ 

r v(*-\)  sinv(x-l)  J 

where  / = - dx 

0 J i xa 

Substitution  of  the  results  obtained  fcr  c,  and  T2  in  Equation  [23, 
gives  the  following  expressions  for  dN/dt  : 


/ ST  + 3\ 

1 3 „ 

3 

3V  + 5 I/2,\  , 

( 2 M 

\2V  + 

2 

8 V’/oK  + aJi 

n/3r  + 3i  f 9 ,3V,  3 V . or  1 
-4»R,B(-j-)e.|J+T+--  VI,-— ij 


t ^bj 


= 3D—(3Y-+-S)\-  •+  — + -V-  - v2!  - ~-l  aoV 
\ <*<  no  2 I in  8 2 '«  4atfj  K,,2 


i 3^c] 


The  integral  1$  is  always  less  than  V " 1 since  t::is  latter  value  's  obtained 
upon  elimination  of  the  x~ 3 sin  F (x  - 1)  factor  of  the  integrand,  this  factor 
being  less  tnan  unity  while  the  remaining  factor  is  always  positive  through- 
out the  range  of  integration.  Moreover  as  seen  in  Figure  ■! , /„  , tr.e  evalu- 
ation of  which  is  discussed  in  Appendix  III,  never  exceeds  C.l  and  is  (2  V ) 1 
for  large  V . Hence  it  is  apparent  that  dN/dt  is  always  positive. 
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EVALUATION  AND  INTERPRETATION  OF  THE  SOLUTION 

The  result  for  JN  t is  illustrated  by  the  curves  shown  in  Figure  1. 
Values  of  l/ty  d77/dt  are  shown  as  a function  cf  K fcr  various  frequencies.  For 
all  of  these  curves  the  values  previously  given  fcr  P0 , D , c„  , n0  ana  pa  were 
used.  Also  used  were  Y = 4/3;pw^l  and  p was  taken  as  1/3.  When  R0  apd 
are  both  extremely  small  then 


v " (£>  »•’)*  - O'  ( j *.f  - [~  B.) 


1 J_ 

2 cm-2 


is  also  small  and  approximately 


1 d X 


~ = 3 D C-^  ( 
dt  un  \ 


3 > + 3 \ / 3 V + 9 


4.52  x 10 


cm2  sec-1  [33] 


Since  for  any  fixed  frequency  «/<^  becomes  small  when  R^  is  suf- 
ficiently small  all  of  the  curves  are  asymptotic  to  this  expression  when  going 
towards  small  R0  . However  it  should  be  noted  that  wnen  the  radius  of  the 
bubble  is  much  smaller  than  10_3c»n,  the  accuracy  of  the  calculation  suffers 
from  the  neglect  of  the  surface  tension  and  for  low  frequencies  from  the  use 
of  the  adiabatic  rather  than  the  isothermal  assumption.  Nevertheless  the  re- 
sults clearly  show  that  extremely  small  bubbles  grew  rapidly  by  diffusion. 

In  fact  the  estimated  rat's  of  growth  of  bubbles  of  small  size  is  so  great  that 
for  low  frequencies  the  growth  per  cycle  is  no  longer  a negligible  fraction  of 
the  initial  bubble  size.* 

♦It  le  clear  that  the  assumption  that  the  growth  of  the  bubble  during  a period  of  the  sound  bees  Is 
negligible  will  Impose  a lower  limit  upon  the  frequenclee  for  which  the  calculation  la  accurate.  The 
aeeuaptlon  as  It  has  been  applied  refers  not  only  to  the  c'lange  In  tae  air  content  duo  to  the  average 
rate  of  growth,  but  also  and  mare  restricting,  to  the  amplitude  of  the  eltematlng  component  of  the  air 
content.  Ir.  particular  It  has  been  assumed  that  this  amplitude  Is  much  smaller  than  the  amplitude  of 
the  volume  variation  l.e.  that  in  amplitude  (N-A/,)/^«I8(R- An  estimate  of  a bound  for  the 
amplitude  of  (N  — AJ)/A^  can  be  obtained  by  considering  the  growth  that  would  occur  In  a half  period  If 
the  conditions  most  favorable  to  growth  l.e.  maximum  expansion  where  constantly  maintained.  Then  at 
the  boundary,  frem  Equation  (18),  talcing  y = 4/8  ; e (/?»„»)  “ where  Rm,,  le  the  maximum 

radius.  Since  steady  conditions  are  assumed  e(r)  “ emURmmJR^T'vl~^tRmm»ir)  + c„ . Hence  from 
Equation  (l6),  dN/dt  m — and  the  roiatlve  growth  in  a half -period,  which 

le  taken  as  an  estimate  of  the  amplitude  of  the  alternat’cg  component  of  the  air  content  Is: 

1 1 /-V  1 _ Zn-A’-,.,  nr c,  21  £«  D_  d 

2 ,\di  v ‘ X 'V  ' Ra  v * *»  n,  !!,?  ft, 

where  A Km,,  * Km,,  ~Ro  and  le  assumed  to  be  small.  Hence  using  the  values  quoted  In  the  text: 


1 1 dS  1 f 3 Alim,,  _ir,  /)  1 6.7  v 10  cm 

2 \ J-  i ■ f Ag  ~ 4 > 0 ft/  v ~ k*v 


If  a ratio  of  l/lOO  is  considered  acceptable,  then  the  aesvmsptlon  Is  satisfied  for  frequencies  such 
that  v > 6.7  jr  10" ‘cm*  . *Biie  le  a rather  conservative  estimate.  It  Is  probable  that  the  as- 

scmcptlon  applies  rulte  well  to  somewhat  lower  frequencies.  Also  It  should  be  noted  that  the  failure 
of  the  assunptlon  fcr  low  frequencies  does  not  Imply  that  the  growth  by  diffusion  would  be  leee  than 
that  calculated;  the  reverse  may  woll  he  the  case. 
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As  R0  increases  towards  resonant  size  the  terms  in  V enter.  Also 
a02  increases.  The  curves  then  depart  from  the  asymptote  on  the  pos’tive  side. 
Naturally  this  occurs  earlier  for  the  higher  frequencies.  As  R0  attains  reso- 
ant  size  m/dt  diverges  to  + «>  because  no  damping  effect  nas  been  included  in 
the  treatment  of  bubble  motion  so  that'  a0  becomes  infinite.  For  bubbles  of 
near-resonant  size  Vis  large,  especially  in  the  case  of  the  lower  frequencies. 
Then  terms  in  V in.  the  bracket  on  the  left-hand-side  of  Equation  [34]  dominate 
Since  for  large  V as  shown  in  Appendix  II  1q  «■  (2F)_1  , an  approximation  for 
1/^  dNjdt  for  the  near-resonant  condition  is 


A _ o n/M  / 3 y +,3\  ( 3 y p„  \t  (R^  /J_\ E d 

Kdt  - 3 Injl  2 I \ 4p„D'M Rj  UW/  R^'  rJ 


(>-&) 


= 0.001815(^2-)^ 


where  Ru  is  the  radius  at  resonance. 

For  the  lower  frequencies  the  damping  is  very  slight  and  it  is  to 
be  expected  that  the  curves  will  in  fact  rise  very  sharply  as  resonance  is 
approached-  However  it  is  difficult  to  estimate  the  height  of  the  peaks  tnat 
will  be  reached  at  resonance-.  According  to  elementary  treatment  of  the  effect 
of  a parallel  sound  beam5  on  an  air-  bubble,  the  rescaant  motion  ei  the  huhblA 
is  a sinusoidal  motion  of  amplitude  pR^/SM,  where  6 's  the  so-called  ^damping 
constant  and  is  a function  of  frequency.  In  the  following  table  the  values  of 
6 are  approximate  empirical  values  read  from  Figure  2,  page  44c>  of  Reference 


Frequency,  cycles  per  sec  . 

1,000 

10,000 

40,000 

Damping  constant,  6 

0.025 

0.11 

0.28 

Nominal  ratio  of  the  amplitude  of  the 

first  harmonic  to  the  equilibrium  radius; 

p/3Yd  for  p = 1/3 

3.3 

0.7 

0.3 

The  relative  amplitudes  show  that  the  elementary  theory  fails  at  low  frequen- 
cies even  for  relatively  small  sound  pressures,  and  the  motion  of  tne  bubble 
is  then  not  sinusoidal  but  very  asymmetrical.*  The  curve's  shown  'n  Figure  1 
have  been  extended  only  to  the  bubble  radii  at  which  ut/u  — 1/2  In  tni s 
case  do  = 1/3  and  for  p — 1/3  the  relative  amplitude  of  the  first  hsrmon’c 

of  the  bubble  motion  is  l/9  so  that  at  tnis  point  the  approximations  made  in 
the  calculation  should  still  be  reasonable.**  The  greater  damping  that  is 

*Under  a*cb  coaditijQ3  one  WOTld  to  find  a nonlinear  amplitude  effect  upon  the  damping 

«When  ufa-  1/2  the  second  harmonic  term  diverges,  i.e,,  a.,  becomes  Infinite.  However  to 
second  order  in  p , a2  does  rot  affect  the  diffusion  and  the  inclusion  of  damping  effects  would 
probably  limit  u,2  to  reasonably  small  values. 


RU  <XTAE  RATE  OF  GRO\  TH  PER  SECOND 
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BUBBLE  RADIUS  ,R0  , CENTIMETERS 

Figure  1 - Relative  Rate  of  Growth  as  a Function  of  Frequency 

and  Bub  ble  Radius 
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observed  for  the  higher  frequencies  would  indicate  that  the  curves  for  frequen- 
cies greater  than  10^  will  not  exhibit  much  of  a peak  st  resonance.  In  foci 
because  of  the  neglect  of  the  damping,  the  rate  of  growth  predicted  for  he 
higher  frequencies  even  in  the  off-resonant  condition  may  be  somewhat  over- 
estimated. It  should  be  noted  that  the  increasing  asymmetry  of  the  motion, 
whether  produced  by  resonant  conditions  or  by  the  sp;.  location  of  sound  beams 
of  greater  intensity,  will  contribute  further  to  tiie  growth  if  tue  bubble  by 
d' ffusion. 

It  is  hard  to  predict  what  happens  when  the  radius  of  tue  bubble 
grows  beyond  resonant  size.  If  the  bubble  becomes  substantially  greater  than 
resonant  size  the  distortion  produced  by  the  motion  may  cause  the  bubble  to 
split  and  thereafter  the  pattern  of  growth  for  small  bubbles  would  aga'n  ap- 
ply. For  very  high  frequencies  the  resonant  bubble  size  is  so  small  that  if 
the  bubbles  do  not  grow  much  beyond  resonant  size  they  would  not  be  evident 
on  visual  inspection.  Also  very  small  bubbles  may  net  have  time  to  rise  cut 
of  the  water  so  that  the  water  would  not  be  deaerated. 

The  time  required  for  a bubble  to  grow  from  one  radius  to  another 
can  be  determined  by  integration  of  the  values  of.  1 /N0  dN/dt  using  tr.e  re- 
lationship dR$/dt  — dN/dt)  . Figure  2 shows  the  result  of  numerical 

integration  applied  to  the  values  for  v = 10  cycles  per  second.  Because  the 
growing  bubole  spends  only  a very  small  part  of  its  life  at  small  radii  it 
would  be  very  difficult  to  determine  tne  original  radius  of  the  bubble  from 
data  on  the  time  after  application  of  a sound  beam  that  is  required  for  a 


TIME  - SECONDS 


Figure  2 - Time  History  of  Bubble  Radius 

This  figure  shows  the  yariation  of  the  radius  of  a bubble  in  a sound  beam  cf  10** 
cycles  per  second.  The  initial  radius  was  taken  as  lo’7  centimeters. 
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APPENDIX  I 

COMPARISON  WITH  BLAKE'S  CALCULATION 

Since  Blake2  includes  the  effect  of  surface  tension  and  uses  an  iso- 
thermal rather  than  an  adiabatic  relationship  between  Pg  and  R his  basic  as- 
sumptions are  more  appropriate  for  verv  small  bubbles  than  the  assumptions 
used  for  the  present  calculation,  the  main  interest  of  which  is  centered  on 
bubbles  of  the  order  of  10~'1  cm  and  greater.  A consequence  of  the  inclusion 
of  surface  tension  is  that  equilibrium  with  respect  to  diffusion  is  not  obtain- 
ed in  the  absence  of  the  sound  beam.  Therefore  Blake  obtains, as  would  be  ex- 
pected, a threshold  sound  pressure  above  which  the  bubble  grows  but  below  which 
the  bubble  dissolves.  However  if  one  is  mainly  concerned  in  obtaining  a real- 
istic value  for  the  threshold  sound  pressure  required  for  the  growth  of  nuclei!, 
it  would  seem  necessary  to  include  in  the  treatment  of  the  problem,  the  influ- 
ence of  whatever  factors  are  responsible  for  the  stability  of  the  nuclei!  with 
respect  to  diffusion.  Unfortunately  definite  information  as  to  the  nature  of 
such  factors  is  not  available. 

It  is  believed  that  the  present  calculation  is  more  detailed  and 
rigorous.  In  the  treatment  of  the  diffusion  problem  Blake  neglects  the  motion 
of  the  boundary  and  also  omits  the  convection  term  in  the  differential  equation. 
Apart  from  the  difference  in  the  motion  of  the  bubble  due  to  the  use  of  the  iso- 
thermal assumption,  it  appears  that  Blake  treats  only  the  off-resonant  case 
since  no  frequency  effect  appears  In  his  equation  for  the  bubble  motion.  Also 
Blake  uses  for  the  most  part  only  first  order  effects.  However  the  ultimate 
result  shows  that  the  growth  of  the  bubble  is  a second  order  effect  in  the 
sound  pressure.  It  is  therefore  advisable  that  second  order  terms  be  retained 
througnout . 

Whereas  Blake  reports  a slight  reduction  of  the  rate  of  growth  of  the 
bubble  when  terms  beyond  the  first  order  in  bubble  motion  are  included,  the  re- 
sults of  tne  present  calculation  show  that  such  terms,  which  reflect  the  a- 
symmetry  of  the  motion,  make  an  Important  contribution.  In  the  off-resonant 
case  where  V is  small,  the  contribution  of  the  second  order  term  as  represented 
by  the  presence  of  aL  ; in  Equation  [J4J,  is  2.5' times  as  great  as  the  contri- 
bution of  the  first  order.  Terms  of  even  higher  order  are  undoubtedly  required 
to  properly  represent  the  effects  of  asymmetry  when  the  resonant  size  is  ap- 
proached. 

When  the  terms  that  reflect  the  lack  of  equilibrium  in  the  absence 
of  the  sound  beam  are  dropped,  and  the  surface  tension  is  taken  as  zero  then 
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APPENDIX  Tl 


THE  EXTRAPOLATION  0?  THE  BOUNDARY  CONDITOK  TC  Rt 

• The  meaning  of  the  expansions  illustrated  by  Equations  [21]  ana  [22] 

may  require  seme  clarification.  Suppose  that  at  time  tj the  radius  of  the 
Dubbie  R is  smaller  than  the  equilibrium  radius  R^  as  shown  in  Figure  3.  Then 
since  the  concentration  of  air  molecules  in  the  water,  car  be  assumed  to  be 
analytic  fo  r > R tr.e  extrapolation  from  R to  R^,  i.e.  Px  to  Pt,  is  clearly 
justified.  However  at  some  later  time  tt  when  the  radius  of  the  bubble  is 
greater  than  R0  .there  is  no  water  at  R0  so  that  physically  the  function,  c 
does  not  exist  at  this  point,  PA  . The  final  justification  for  these  expansions 
is  simply  that  by  applying  them  a solution,  e(r,  t ) , is  obtained  that  is  ana- 
lytic for  all  values  of  r> 0.  It  is  certainly  permissible  to  apply  this  so- 
lution for  values  cfr^Rand  it  ’s  apparent  that  the  setting  of  the  boundary 
conditions  at  RQ  in  the  manner  described  raaKes  certain  that  the  solution 
will  satisfy  the  proper  boundary  conditions  at  R. 


c of  t. 


Figure  3 - The  Extrapolation  of  Boundary 
Conditions  from  R to  R^ 
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APPENDIX  III 
EVALUATION  OP  /« 

The  Integral  ^(V)  = J S-n  V ^ — -dz  can  be  expressed  in  t»erms 

of  the  generalized  expohential  Integrals.  Integrating  by  parts  so  as  to 
reduce  the  exponent  of  the  denominator  one  obtains: 

. . \ 

I(V)  =|  - ev  V|[J  -EAO.V)]  sin  V +[Ee  (0 ,V)  -Y-  InV-  ^ in2jcosVj 

where  V ». 5772157  ...  Is  Euler's  constant  and  EJO,  V) , £^(0,  V)  are  the  gener- 
alized exponential  Integral  functions.  These  Integrals  have  been  tabulated®. 
For  values  of  V smaller  than  the  tabulated  range  the  following  series 
expansions  can  be  applied. 

Vs  V4 

Ec(0,V)  = V-  J+24  + -'- 
V2  V3 

EA0,V)  * v - Y + y 

For  values  of  V greater  than  the  tabulated  range  the  integral  70  can  be  evalu- 
ated directly  using  the  expansion 

which  can  be  obtained  either  by  Integrating  by  parts  so  as  to  increase  the 
exponent  of  the  denominator  of  the  integrand,  or  by  expanding  the  denominator 
Vx*  in  a power  series  about  s * 1.  For  values  of  V that  are  Integral  multi- 
ples offfjthe  integral  can  be  conveniently  evaluated  by  use  of  the  rapidly 
converging  series 

I0(mn)  = V2  T*  Kl(n  + m)n]e~"w 

AWn  = 0 


A plot  of  /0 ( V ) over  a considerable  range  of  the  variable  is  shewn  in  rigure 

4. 


\ 


V 


Figure  4 - Plot  of  4<V)  = J, 


i j sin  V (z  - 1)  ^ 

z' 


1 


r € ~ 8 gjp  x 

where  values  of  K(v)  = Jo  ~ ~ dz  are  readily  obtained  by  numerical  lnte 
gratlon.  In  this  manner  the  following  values  were  obtained. 


V 

/. 

V 

/. 

n 

0.07393 

ln 

0.01991 

2 n 

0.09186 

0 n 

0.01771 

3 n 

0.03944 

9n 

0.01994 

4 TT 

0. 03172 

10  n 

0.01449 

5 n 

0.02649 

11  n 

0.01328 

6;'i 

0.02274 

12  n 

0.01227 
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